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Abstract 

The effect of the position of on-ramp and off-ramp i\ and %i, respectively, on the 
one dimensional-cellular automaton traffic flow behaviour, is investigated numerically. 
The injection rates at i\ and i% are ao an d (3o, respectively. However, in the open 
boundary conditions, with injecting and extracting rates a and (3 and using parallel 
dynamics, several phases occur; namely, low density phase (LDP), intermediate density 
phase (IDP), plateau current phase (PCP) and high density phase (HDP). It is found 
that phase diagrams exhibit different kind of topologies. For intermediate value of ex- 
tracting rates (3q and (3 and low value of a, (ii,ao) phase diagram provides LDP-IDP, 
LDP-PCP, IDP-PCP, and PCP-HDP transitions, and critical end points. The off-ramp 
position is located to the middle of the road. By increasing (3q and /?, the IDP desap- 
pears. For high value of (3, only LDP-HDP persists. 

Pacs number : 05.40.-a, 05.50. +q, 64.60. Cn, 82.20.wt 

Keywords : Traffic flow, phase diagrams, Cellular automata, on- and off-ramp, parallel 
update. 

*corresponding author e-mail address: mhirech@fsr.ac.ma 



1 



1 Introduction 



During the last years, the field of transport have attracted several researchers [ ]. This 
interest is due primarily to the fact that transportation problems are related to the 
global behaviour of systems with many elements interacting at short distances, such as 
the vehicles traveling on the streets, or informations which travel over the internet net- 
work. In particular, the investigation of open traffic systems with on- and off-ramps is 
quite popular at the moment [ ] . One reason for this is the impact of the understanding 
of varying the different flow ratesin order to optimize the total flow or trip times. 

Among the different methods of investigation and simulation of highway traffic, as- 
symetric simple exclusion process (ASEP) is the most promesing [ ]. Indeed, ASEP is 
the simplest driven diffusive system where particles on a one-dimensional lattice hop 
with asymetric rates under excluded volume constraints. 

The question we want to answer is the following: let us suppose we have a highway 
running in an urban conglomeration and that there are an access from urban conglom- 
eration to the highway and an exit from the highway. We want to understand where 
the access and the exit positions must be located in order to maximize the flux of cars 
in the road. Our aim in this paper is to study the effect of the on-ramp and off-ramp 
positions on the one dimensional-cellular automaton traffic flow behaviour in the open 
boundaries case. Depending on the injecting and extracting rate values, an adequate 
localization of the on- and off-ramp positions leads to the appearance of new phases 
and topologies. Moreover, to compare our results to those where only one off-ramp was 
taken into account [30], quantitative differences can be understood from the behaviour 
of average density, current and phase diagrams for different parameters. 

The paper is organised as follows: Model and method are given in section 2; section 
3 is reserved to results and discussion; the conclusion is presented in section 4. 
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2 Model 



We consider a one-dimensional lattice of length L. Each lattice site is either empty or 
occupied by one particle. Hence the state of the system is defined by a set of occupation 
numbers ti,T2,...,tl, while Tj = 1 (ti = 0) means that the site i is occupied (empty). We 
suppose that the main road is single lane, an on-ramp and an off-ramp connect the main 
road only on single lattice i\ for entry and on single lattice 12 for way out. During each 
time interval At, each particle jump to the empty adjacent site on its right and does not 
move otherwise [i 7^ 12)- At is an interesting parameter that enables the possibility to 
interpolate between the cases of fully parallel (At = 1) and random sequential (At — > 0) 
updates [29]. Particles are injected, by a rate a At, in the first site being to the left 
side of the road if this site is empty, and particles enter in the road by site i±, with a 
probability aoAt without constraint, if this site is empty. While, the particle being in 
the last site on the right can leave the road with a rate (3 At and particles removed on 
the way out with a rate o At. At site i\ (12) the occupation (absorption) priority is 
given to the particle which enter in the road (particle leaving the road). Hence the cars, 
which are added to the road, avoid any collision. 

In our numerical calculations, the rule described above is updated in parallel, At = 1, 
i.e. during one update step the new particle position do not influence the rest and only 
the previous positions have to be taken into account. During each of the time steps, 
each particle moves one site unless the adjacent site on its right is occupied by another 
particle. The advantage of parallel update, with respect to sublattice or sequential 
update is that all sites are equivalent, which should be the case in realistic model with 
translational invariance. 

In order to compute the average of any parameter w (< w >), the values of w(t) obtained 
from 5 x 10 to 10 5 time steps are averaged. Starting the simulations from random 
configurations, the system reaches a stationary state after a sufficiently large number 
of time steps. In all our simulations, we averaged over 60 — 100 initial configurations. 
For the update step, we consider two sub steps as shown in figure 1: 
In the first sub step, the sites i\ and 12 are updated and in the second half, the chain 
updates. Thus if the system has the configuration T\(t), T2(t),...,TL(t) at time t it will 
change at time t + At to the following: 
For i = i\ , 

n(t + At/2) = i (l) 

with probability 

q t = n(t) + Ml - n(t)) - n(t)(l - n +1 (t))]At (2) 

and 

T;(t + At/2) = (3) 
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with probability 1 — ft. Where i\ and «o denote the position of the entry site and the 
injection rate, respectively. 
For i = i 2 , 

n (t + At /2) = 1 (4) 

with probability 

ft = Tj(t) + - Ti) - #)Ti(t)]At (5) 

and 

n (t + At /2) = (6) 

with probability 1 — ft. Where ^2 and /?o denote the position of the absorbing site and 
the absorbing rate, respectively. 
For l<i<L with i ^ i\ and i ^ i 2 , 

Ti(t + At) = l (7) 

with probability 

ft = n(t) + [7i_i(t)(l - Ti(t)) - Ti(t)(l - r i+ i(t))]A/ (8) 

and 

n(t + At) = o (9) 

with probability 1 — ft. 
For i=l, 

Ti(t + At) = l (10) 

with probability 

ft = n(i) + [q(1 - n(t)) - n(t)(l - T 2 (t))]At (11) 

and 

n(t + At) = (12) 

with probability 1 — ft. 
For i = L, 

r L {t + At) = l (13) 

with probability 

QL = r L (t) + [r L _i(t)(l - r L (i)) - Pr L (t)]At, (14) 

and 

r L (t + At) = (15) 

with probability 1 — ft,. 
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3 Results and Discussion 



As we have mentioned previously, our aim in this paper is to study the effect of the 
positions of the on- and off-ramps i\ and 12, respectively, for different values of a, (3q and 
(3, on the average density and flux in chain. The study is made in the open boundary 
conditions case, ao and a denote the injecting rates at first site (i = 1) and at site i±, 
respectively. Pq and (3 are the extracting rates at site 12 and at the last one (i = L), 
respectively. The length of the road studied here is L=1000. 

The figures 2(a) and 2(b) give respectively the variation of the average density p and 
average current J versus the injecting rate ao for several values of the on-ramp position 
i\. These figures are given for a = 0.1, (3 = 0.1 and (3q = 0.4. However, when the 
on-ramp is located upstream of the off-ramp, the system studied exhibits four phases, 
depending of the behaviours of the density, p, and the current, J. Namely: i) The low 
density phase (LDP), where the averages density and current increase when increasing 
the rate of injected particles ao- h) The intermediate density phase (IDP) characterised 
by a smoothly increase of the density and average current, iii) The plateau current 
phase (PCP) for which the density and current are constant in a special interval of 
ao- iv) The high density phase (HDP) in which, for high values of ao, the current de- 
creases and the density reaches its maximum value and remains constant. On the other 
hand, when i\ is located downstream from 12, the IDP desappears. In addition, when 
increansing i\ for a given value of ao, the figure 2(a) shows that the average density 
is constant in LDP, increases in IDP and PCP then decreases in HDP. The figure 2(a) 
exhibits an inversion point situated at the PCP-HDP transition. Moreover, the figure 
2(b) shows that the average current deceases by increasing i±, for any value of ao- 
For i\ < i2, we note that the IDP, which doesn't appear in the model where only the 
off-ramp is taken into account [30], occurs for the intermediate values of ao (ao c i < 
ao < aoc2)- «0ci and ao C 2 correspond to the transition between LDP-IDP and IDP- 
PCP, respectively. While the PCP arises between two critical values ao C 2 and ao C 3 of 
injecting rate ao- Where ao C 3 corresponds to the PCP-HDP transition. Note that the 
transitions which occur at ao C 3 disappears when i\ is located after 12 (Figure 2a). Now, 
in order to have a suitable criterion for determination of the nature of the transition, 
we identify the first order transition (abrupt transition) by the jump in the average 
density or by the existence of a peak in the derivative of p(ao) with respect to ao- The 
jump in density corresponds to a first order transition [29]. This means that the above 
transitions are of first order type. 

Collecting the results illustrated in figures 2(a) and 2(b), the four regions are given 
on the phase diagram (ii,ao) shown in Figure 2(c). Beside this, such phase diagram 
exhibits four critial end points, around which there is no distnction between the phases. 
This critical end points are indicated by CEP. 
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For low values of a, f3 and (3o (a = 0.1, j3 = 0.1 and fto = 0.1), The (i\,ao) phase dia- 
gram is presented in figure 3. This figure exhibits only tow first order phase transitions. 
Namely, LDP-PCP, PCP-HDP transitions. The later one can be founed by variying 
ao^or a given value of i\ lower than 12, or at i\ = 12, for 0.45 < uq < 0.95. Moreover, 
the figure 3 exhibits two critical end points. The comparison of figures 2 and 3 highlight 
the effect of (5q on the ao) phase diagram. Indeed, for intermediate value of (3$, the 
IDP arises. 

For a sufficiently large value of f3o, the critical end points desapears, as shown in figure 
4. This figure is given for a = 0.1, (3 = 0.1 and [3q = 0.8. 

The figure 5 give the (ii,ao) phase diagram for a = 0.1, (3$ = 0.1 and (3 = 0.3. In 
this case, the system exhibits four phases and three critical end points. From figures 2c 
and 5, we deduce that the IDP arises for intemediate values of extracting rates (3 or (3q, 
when i\ is upstream of 2- 



4 Conclusion 



Using numerical simulations, we have studied the effect of the on- and off-ramp posi- 
tions on the traffic flow behaviour of a one dimensional-cellular automaton, with parallel 
update. Depending on the values of a, (3 and fio, the («i,ao) phase diagram exhibits 
different topologies. The IDP occurs only at special positions of on- and off-ramps with 
special values of extracting rates (3 and (3q and injecting rates a and ao. The transition 
between different phases are of first order. Furthermore, the system exhibits critical 
end points in the (ii,oco) plane in the case of moderate values of and (3q and small 
value of a. 
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Figure captions: 

Fig.l: Example of configuration obtained after two sub steps for system size L = 15. 
Fig.2: For a = 0.1 and (3 = 0.1; (a) Average density p versus the injection rate oq; (b) 
Variation of the average current as a function of Qo; (c) Phase diagram (/?o,ao)- The 
number accompanying each curve, in (a) and (b), denotes the values of (3q. 
Fig.3: Phase diagram (/3,ao) for a = 0.1 and (3q = 0.1. 

Fig.4: Phase diagrams (/3o,«o) for (3 = 0.1; (a) a = 0.2; (b) a = 0.4; (c) a = 0.5. 
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